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ABSTRACT 

This paper describes the flywheel test facility developed at 
the NASA Glenn Research Center with particular emphasis on 
the motor drive components and control. A 4-pole permanent 
magnet synchronous machine, suspended on magnetic 
bearings, is controlled with a field orientation algorithm. A 
discussion of the estimation of the rotor position and speed 
from a "once around signal" is given. The elimination of small 
dc currents by using a concurrent stationary frame current 
regulator is discussed and demonstrated. Initial experimental 
results are presented showing the successful operation and 
control of the unit at speeds up to 20,000 rpm. 

INTRODUCTION 

One of the key components of the flywheel energy storage 
system is the electric motor and its control. Energy storage and 
recovery are achieved by using the motor to increase or 
decrease the flywheel rotor speed as necessary. Good control of 
the motor is thus very important for the proper operation of the 
flywheel system. As part of the flywheel technology 
development effort, NASA Glenn Research Center has built a 
test facility with the capability to rapidly test and evaluate 
advanced motor control algorithms. It is the purpose of this 
paper to describe the test facility (with particular emphasis on 
the motor control portion), the basic motor control algorithms 
developed and to present initial experimental results. 

NOMENCLATURE 

Lq is the q-axis machine inductance, henries, 
is the d-axis machine inductance, henries, 
q-axis voltage or current in the rotor reference frame. 

fj d-axis voltage or current in the rotor reference frame. 

f* q-axis voltage or current in the stator reference frame. 


d-axis voltage or current in the stator reference frame. 

p is the derivative operator, d/dt. 

0 r is the angle between the stator q-axis and the rotor q- 
axis, radians. 

X, a f is the flux linkage due to the rotor magnets, volt-sec. 
is the electrical rotor speed, radians/second. 

FLYWHEEL TEST FACILITY 

The flywheel system test facility consists of a test cell and 
a control room. The test cell is physically separated from the 
control room for safety purposes. The present flywheel system 
under test consists of a 4 pole permanent magnet synchronous 
machine, a high strength composite rotor, magnetic bearings 
and a housing structure which is sealed and pumped to a low 
vacuum with water cooling capability. The flywheel system 
itself is housed in a containment structure that is closed during 
operation. Power for the motor is derived from a standard six 
switch three phase inverter with a dc power supply source. The 
inverter is a commercial off the shelf intelligent power module 
(gate drive circuitry included in the unit) rated at 600 volts, 200 
amps and uses IGBTs for the power switches. 



Figure 1 : Block Diagram of Flywheel Motor Control 
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The control room has independent stations for various 
functions: magnetic bearing control, motor control dc power 
supply control and temperature monitoring. The motor control 
function is implemented digitally using a commercially 
available microprocessor board. The control algorithms are 
written in a high level block diagram simulation language and 
converted to the necessary microprocessor code via software 
provided by the manufacturer of the microprocessor board. 
The motor current is measured in the test cell and fed back to 
an analog to digital (A/D) converter board in the control room. 
The inverter power switch duty cycles are calculated digitally 
as part of the control algorithm and output as analog signals 
using a digital to analog (D/A) converter. The duty cycle 
analog signals are sent to a pulse width modulation (PWM) 
generation board in the test cell. The output of the PWM 
generation board is a set of gate drive signals for the inverter 
unit. A block diagram of the motor control is shown in Figure 
1 and described extensively in the next section. 

MOTOR CONTROL ALGORITHM 

The motor control algorithm is based on the principles of 
field orientation control or ’’vector control” [1]. The basic idea 
behind this control technique is to orient the applied stator 
currents to the rotor magnetic field. When this is done, the 
motor control is simplified because the control variables 
become dc quantities in steady state. The operating point of the 
machine can then be accurately and dynamically controlled so 
that high efficiency and fast response are obtained. A brief 
mathematical description for this technique follows. 

A three phase machine, without a neutral connection, can 
be equivalently described as a two phase machine through a 
transformation from abc coordinates to dq coordinates as 
follows [2]. 


II 

(1) 

1 2 

fd = V3 fa 'V3 fb 

(2) 

The reverse transformation is 

fa=fq 

(3) 

f _ I , y[l f 

“ “2 m ‘ 2 

(4) 

1 x/3 

fc = -2 f q + 2 f d 

(5) 


The d and q variables described in (1) and (2) are in the 
stator reference frame. From (1) it can be seen that the q-axis 
is aligned with the ’a' phase. This means that the 'a’ phase 
current is equal to the q-axis current and the ’a' phase voltage 
(Van) is equal to the q-axis voltage. 

In a permanent magnetic machine, it is convenient to 
transform these variables to a reference frame that is rotating 
synchronously with the rotor magnetic field. In the rotor 
reference frame, the d-axis is defined to be co-linear with the 
rotor magnetic field axis. The transformation from the stator 
frame to the rotor frame is given by (6) and the inverse 
transformation is given by (7) [2]. 



The permanent magnetic synchronous machine can then be 
modeled in the rotor reference frame as follows [1]* Equation 
(9) gives the stator flux linkages in the rotor reference frame 
and (10) is the torque expression. 


^qs _ Rs + LqP oyL^ Tijs r 

V<j s _ _ -0\-Lq Rs + EdpX*ds_ L 0 J 


( 8 ) 



4 [\is ^qs " ^qs *ds] 


In field orientation control the d-axis current, *dV 1S 
commanded to 0. From (9) and (10) it can be seen that this 
results in a simplified expression for torque, given in (1 1). 



( 11 ) 


Because the magnitude of X a f is constant due to the rotor 
magnets, the torque of the machine is proportional to the q-axis 
current. This result is similar to the dc motor with a constant 
field winding current where the torque is equal to the armature 
current times the torque constant. Thus control of the torque is 
achieved by properly controlling the rotor reference frame 
currents, i q r s and ij s . 

The torque command to the machine is determined from 
the outer loop control. In the results presented here, the outer 
loop is a speed controller, that is, errors in speed result in a 
correcting torque command. In an actual flywheel energy 
storage system, the torque command is derived from dc bus 
current or voltage commands. This is described in a 
companion paper [3]. 

To control the currents in the machine to be the 
commanded values, a current regulated voltage source inverter 
is used. This means that current errors result in voltage 
commands to the inverter which increase or decrease the 
applied voltage to increase or decrease the current respectively. 
This is accomplished through a synchronous frame current 
regulator [4] that is basically a PI controller operating on the 
rotor reference frame currents. The output of the controller is 
considered to be a voltage command in the rotor reference 
frame that is then transformed to the stator reference frame. 
The stator frame voltage commands are then used to calculate 
the inverter switch duty cycles. 

There are several methods to find the duty cycles from the 
commanded stator frame voltages [5]. In the implementation 
used here, space vector modulation is used. Space vector 
modulation is a digital technique to calculate the duty cycles 
directly from the stator reference frame d and q voltages. One 
advantage of space vector modulation is that it increases the dc 
bus utilization to the maximum value. This means that for a 
given dc bus voltage, the maximum fundamental phase voltage 
is achievable by using space vector modulation. This is in 
contrast to the more common sine-triangle modulation which 
results in a lower dc bus utilization. 

In this implementation, the duty cycles are output as 
analog signals which are sent to a PWM generation board. The 
analog signals are compared to a carrier signal triangle wave to 
create the gate drive signals to the inverter. This method allows 
the switching frequency to be separate from the sample 
frequency which is limited by the calculation time of the 
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control algorithm. In this implementation, the sample time is 
set to 100 usee and the PWM switching frequency is set to 
twice that value at 20 kHz which is the maximum switching 
frequency of the inverter. 

It can be seen by the rotor reference frame transformations 
of (6) and (7) that the rotor reference frame quantities are based 
on the rotor angle, 0 r . Thus it is necessary' to know what this 
angle is. In industrial implementations of field orientation, it is 
common to use an encoder or resolver to feed back the rotor 
position. The high speed nature of this application makes it 
difficult to use the traditional feedback devices. Instead, the 
rotor position is estimated based on knowledge of a "once 
around" (OAR) position marker. In this machine, one half of 
the rotor hub is coated with a dark oxide, the other half is 
polished aluminum. The once around signal is based on the 
optical detection of the dark or light surface and thus a square 
wave at the same frequency as the machine speed is generated 
with the rising edge corresponding to the change from dark 
oxide to polished aluminum. Due to the large inertia of the 
flywheel rotor, the speed is reasonably assumed to be constant 
between once around signals and from this, the position can be 
estimated. This will be discussed further in the next section. 

INITIAL IMPLEMENTATION ISSUES 

Back EMF Test 

As explained previously, to properly control the motor 
using the field orientation technique, the instantaneous position 
of the rotor q-axis with respect to the stator q-axis (0 r ) must be 
known. The once around signal indicates the position of the 
color change (dark to light) on the rotor but this does not 
necessarily coincide with the position of the rotor q-axis. What 
needs to be determined is the angle between the rotor q-axis 
and the stator q-axis as a function of the once around signal. 
This can be determined experimentally by measuring both the 
phase voltages of the machine and the once around signal as the 
rotor freely decelerates ("spins down"). For a freely 
decelerating machine, the stator current is 0 so the phase 
voltages are entirely due to the back emf of the machine. 
Stated mathematically, (for 0 stator current) [1]: 

V qs = P^qs where \f s = V in < 0 r) ( 12 > 

V ! = P^k where ( |3 > 

By definition, because the q-axis in the 

transformation of (1) is aligned with the stator a-phase. Thus in 
a freely decelerating machine, the relationship between 0 r and 
Van is given in (14). 

=^08(0,) (14) 

V^ can b e measured if there is a neutral terminal on the 
machine or V ab can be measured and V^ calculated. In this 
case, the motor neutral was available for measurement purposes 
(not normally connected). Figure 2 shows a plot of the once 
around signal and the V^ voltage. 

From (14) it can be seen that when 0 r = 0, V an will be at a 

maximum. Figure 2 shows that this maximum occurs twice in 
one cycle of the OAR signal at approximately the middle of 
each pulse. The machine is a 4 pole machine so between every 
rising edge of the OAR signal, 360 mechanical and 720 


electrical degrees will have passed. Figure 2 shows that the 
V an maximum occurs approximately 90 mechanical degrees or 
180 electrical degrees behind the rising edge of the OAR 
signal. This means that the 0 r transformation angle is equal to 
the electrical angle determined bv the OAR signal minus 180°. 



Figure 2: Phase voltage and OAR signal for freely decelerating rotor. 

Once the relationship between 0 r and the once around 
signal is known, the transformation from the stator to the rotor 
reference frame can be done. From (8), it can be seen that for a 
freely decelerating rotor, when i q r s = i d r s = 0, the back emf 
voltage appears entirely in the rotor reference frame q-axis and 
V d r s = 0. This means that the transformation angle 0 r can be 
fine-tuned. The voltages measured in the back emf test can be 
transformed to the rotor reference frame using the relationship 
0 r = OAR electrical angle - 180°. (15) 

If the transformation angle 0 r given by (15) is not exactly 
correct, then VJ S will not be exactly zero. Figure 3 shows a 
plot of V d r s where 0 r lags the OAR electrical angle by 1 80°. It 
can be seen that V d r s is slightly less than zero. Figure 3 also 
shows the case where 0 r lags the OAR electrical angle by 170°. 
It can be seen that this average value of V d r s is now more nearly 
equal to zero. Thus the actual transformation angle used in the 
control algorithm is given by (16). 

0 r = OAR electrical angle- 170° degrees (16) 

The ac components of the q and d axis voltages in Fig. 3 
are due to the spatial harmonics of the motor windings and are 
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Figure 3: V^ for Variations in 0 r 
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Starti ng the machine 

To operate the machine using field orientation, the rotor 
position angle must be known. However, the machine must be 
spinning to estimate the position from the once around position. 
Initially, this dilemma was resolved by starting the machine in 
"open loop" operation. A fixed voltage and frequency were 
applied to the machine directly until the rotor started to move. 
Once a OAR signal was established, the control of the machine 
was shifted to closed loop speed control. 

It was additionally observed that the machine could start 
under closed loop control if a low speed command was given. 

EXPERIMENTAL RESULTS 

Closed Sp eed Loop Operation 

Closed loop speed control was achieved by using a PI loop 
around the torque control as shown in Fig. 1 and the unit was 
tested at speeds up to 20,000 rpm. Using the once around 
signal, two methods were developed to estimate the speed and 
position. The first method works w>ell at low speeds and 
degrades as the speed increases due to quantization errors. The 
second method is limited at low' speeds due to saturation of the 
digital counter register but works very well at high speeds. 

Speed Estimate: Method One 

In this method, the speed was determined by counting the 
number of sample periods, t p , between once around signals and 
multiplying by the sample rate, T s as shown in (17). The speed 
was assumed to be constant between OAR signals leading to 
the position estimate shown in (18) where t pz is the number of 
sample periods since the last rising edge of the OAR signal. 
The position estimate was reset to zero at every rising edge of 
the OAR signal. 


„ 60 
Speedjp m - 

(17) 

2 k 

® electrical 1 ^ = ^P^^rpm t pz 

(18) 


The sample rate of the control algorithm is 100 usee. This 
leads to the possibility of error shown in Table 1. The 
quantization errors of this method led to speed oscillations 
which became progressively worse as the speed increased. 
Figure 4 shows the speed estimate for a commanded speed of 
5000 rpm. Oscillations on the order of 12 rpm are clearly 
evident. Figure 5 shows the impact of the speed oscillations on 
the phase current. 


Actual 
Machine 
Speed (rpm) 

Maximum Possible 
Speed Estimate Error, 
Method One, rpm 

Maximum Possible 
Speed Estimate 
Error, Method Two, 
rpm 

1,000 

1.66 

.03 

5,000 

41.32 

.83 

10,000 

163.93 

3.33 

20,000 

645.16 

13.32 

60,000 

5454.54 

119.76 


Table 1 : Possible Speed Errors for Alternate Methods 


Speed Estimate: Method Two 

This method w'as based on a commercially available timer 
board which contains two 16 bit counters. Counter 1 is used to 


5010 
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4995 


4990 

0 0*2 0.4 0.6 08 1 

time, seconds 

Figure 4: Speed Oscillation due to Quantization Error with Method 
One Speed Estimate 


2.5 



Figure 5: Oscillations in Phase Current due to Speed Estimate 

Oscillations for Method Two 

determine the OAR frequency and counter 2 is used to 
determine the time that has elapsed since the most recent rising 
edge of the OAR signal. The counter clocks are scaled so that a 
signal in the range of 7.64 Hz to 500 kHz may be detected. 
(The clock rates may be scaled higher in software as necessary 
for higher rotor speeds.) The speed is determined from counter 
1 and (19) and the OAR electrical angle is determined from 
(20) where T is the output of counter 1 and Y is the output of 
counter 2. 

f 

Speedjpm = 60 2 (19) 

® electrical = 4rcft (20) 

The motor speed using speed estimate method two for 
feedback is shown in Fig. 6. The oscillations in both the speed 
estimate and the phase current (Fig. 7) are eliminated. 

DC Offset Currents 

Figure 8 shows the frequency spectrum of the phase 
current at an operating speed of 10,000 rpm. The fundamental 
component, at 333 Hz, clearly dominates. However, there is a 
significant dc component also present. Although not predicted 
theoretically, the dc component may be caused by slight 
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Figure 6: Speed oscillation with Method Two Speed Estimate 


operates in the stationary reference frame. This current 
regulator is a PI controller tuned to have a very low bandwidth 
so it essentially only regulates the dc current as shown in Fig. 9. 
Due to the low bandwidth of the PI, the fundamental current 
passes through the regulator unchanged while the dc current is 
regulated to zero. Using the stationary frame current regulator 
in addition to the synchronous frame current regulator resulted 
in the elimination of the dc component as shown in Fig. 10. 



Figure 9: Block Diagrams of d- and q-axis Stationary Frame DC 
Current Regulators 
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Figure 7: Phase Current with Method Two Speed Estimate 



Frequency, Hz 

Figure 8: Phase Current Frequency Spectrum at 10,000 rpm 

mismatches or offsets in the phase voltages applied to the 
machine. 

The synchronous frame current regulator was tuned for a 
bandwidth of 200 Hz. The synchronous electrical speed at 
10,000 rpm is 333 Hz. This means that currents in the 
frequency range of 133 Hz to 533 Hz are within the bandwidth 
of the current regulator. Thus the current regulator will not 
effectively eliminate the dc currents. One approach to this 
problem is to introduce an additional current regulator which 


0.8 -r T v 



Frequency, Hz 


Figure 10: Phase Current Frequency Spectrum at 10,000 rpm with 
stationary frame current regulator 

Sensitivity to Error In Position Estimate 

In this implementation of field orientation control the 
commanded d-axis current, i d r s , is set equal to zero. In this 
strategy, the direct axis stator flux linkage, is constant and 
due entirely to the rotor magnet with no contribution from the 
stator current. The torque is produced by the interaction of this 
flux with the q-axis current according to (11). However, if the 
rotor reference frame transformation angle, 0 r , is not correct, 
then there will be d-axis current in addition to the commanded 
q-axis current. For small errors in angle estimate, this 
additional current does not contribute significantly to the total 
phase current in the machine. However, as the error 

approaches 90 electrical degrees, the theoretical maximum, the 
total stator current increases rapidly. This is shown in Fig. 11. 
The data points in Fig. 1 1 were taken at a speed of 7,000 rpm. 

20000 RPM OPERATION 

The speed of the motor was limited to 20,000 rpm due to 
limits on the particular flywheel rotor used in the system. The 
following results show steady state operation at 20,000 rpm. 
The small magnitude of the phase current (rms value just over 1 
amp) indicates the very low losses of the machine when 
suspended on magnetic bearings and operated in a vacuum. 
Fig. 12 shows the motor phase current. The current looks 


NASA/TM — 2001-210986 


5 


Copyright © 2001 by ASME 




somewhat distorted however an examination of the spectrum 
shown in Fig. 13 shows the major frequency component at 667 
Hz as expected. Also, the current shown in Fig. 12 has an 
extremely small magnitude which may contribute to the 
perceived distortion. 


8 - — 



°0 20 40 60 80 100 

Number of Degrees of Rotor Angle Error 

Figure 1 1: Sensitivity of Phase Current to Rotor Angle Error 



Figure 12: Motor Phase Current at 20,000 rpm Operation 



Figure 13: Motor Phase Current Spectra for 20,000 rpm Operation 


Finally, Fig. 14 shows the duty cycle for the phase A high 
switch at 20,000 rpm with a 120 volt dc bus for both sine- 
triangle modulation and space vector modulation. This plot 
shows that the space vector modulation technique needs less dc 
bus voltage than the sine-triangle technique to generate the 
same phase voltage. The space vector duty cycle is seen to be 
well within the duty cycle range of zero to one for operation at 
this speed from this voltage. However, the duty cycle for the 
sine-triangle technique is clamped at one at some points and at 
zero at others. This duty cycle which would result in a higher 
distortion of the synthesized phase voltage. 



Figure 14: Phase A High Switch Duty Cycle for Space Vector 

Modulation and for Sine -Triangle Modulation 

CONCLUSIONS 

This paper has presented the initial experimental results of 
the flywheel motor test facility at the NASA Glenn Research 
Center. The once around signal was successfully used for 
speed and position feedback. The back emf test and the 
subsequent fine-tuning was shown to accurately determine the 
relationship between the OAR and the rotor q-axis. The dc 
offset current was eliminated with the addition of the 
concurrent stationary' frame regulator. Space vector 
modulation was shown to increase the dc bus utilization over 
sine-triangle modulation. 
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